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Abstract 

Since the spin of real particles is of order of h, it is difficult to distinguish in a 
quantum mechanical experiment involving spinning particles what part of the 
outcome is related to the spin contribution and what part is a pure quantum 
mechanical effect. We analyze in detail a classical model of a nonrelativistic 
spinning particle under the action of a potential barrier and compute numeri¬ 
cally the crossing for different potentials. In this way it is shown that because 
of the spin structure there is a nonvanishing contribution to crossing for en¬ 
ergies above a certain minimum value, even below the top of the potential 
barrier. Results are compared with the quantum tunnel effect. 
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I. INTRODUCTION 


The aim of this work is to describe in classical terms the classical spin contribution to 
the so called tunnel effect by assuming a particular spin structure of a classical spinning 
particle. 

In section we make some general considerations about the classical limit of quantum 
mechanics, and also about the definition of elementary particle in quantum mechanics and 
the different treatment deserved to spin in classical as well as in quantum mechanics. These 
considerations lead us to propose a general definition of elementary particle at the classical 
level that resembles the one used in quantum mechanics. 

A Lagrangian approach of a nonrelativistic spinning particle model to be used later, 
is described in detail in Sec. 0 where the main features are that the particle spin arises 
because the separation between its center of mass and center of charge. The center of mass 
satisfies Newton’s equations where the external force is defined on the point charge position 
while this center of charge moves around the center of mass with an isotropic harmonic 
motion giving rise to the magnetic properties of the particle. The canonical approach of this 
model is also presented in Sec. B 

We shall consider the solutions of a one-dimensional model of this classical spinning 
particle interacting with a potential barrier in Section |V|. Dynamical equations are nonlinear 
and a general solution is dificult to obtain in closed form even when the external potential 
is a linear function of the charge position. We solve the problem by numerical integration 
showing the penetration effect for different initial conditions and different slopes of the 
potential barrier. The crossing is achieved because of the separation between the center 
of mass and center of charge in such a way that the average repulsive force acting on the 
spinning particle is smaller than in the point particle case. Results are compared with the 
quantum tunnel effect for a point particle under the same one-dimensional potential. 
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II. CLASSICAL LIMIT OF QUANTUM MECHANICS 


One usually reads in textbooks on quantum mechanics that Classical Mechanics (CM) 
is the limit h —>■ 0 of Quantum Mechanics (QM) |^. Of course there are exceptions. For 
instance in Feynman’s Lectures on Physics [Q we read: ‘In the classical limit, the quantum 
mechanics will agree with Newtonian Mechanics’. Both expressions are equivalent if by 
CM we understand the classical mechanics of spinless (or point) particles, i.e., Newtonian 
Mechanics. We think this is the meaning the mentioned authors try to express. But the hrst 
statement as it stands, might lead to wrong interpretations if considered literally, as some 
of the quoted references might suggest. In order to clarify this idea of the classical limit of 
QM, let us consider the following simplihed diagram of Fig. Q. 

The set A represents the whole body of knowledge of CM which includes two subsets, 
the subset (s = 0) or domain of spinless particles or Newtonian Mechanics and the subset 
(s 7 ^ 0) of spinning particles. If we restrict CM to satisfy the additional requirement of 
the Uncertainty Principle we enter into the more restricted body of knowledge of QM, 
represented by the smaller set B in which we also have the two subdomains of spinless 
and spinning particles. If in this more constrained domain we perform now the additional 
restriction of taking the limit h ^ 0 it is out of logic that after these two restrictions we 
shall reach the wider and less restrictive domain of the whole CM. 

In QM we have that the measurement of any two observables C and D is not in general 
compatible, and the uncertainty in their simultaneous measurement is related to its commu¬ 
tator [C, D\ cs h wich is of order of Planck’s constant h. But we also have in QM that the 
particle states satisfy eigenvalue equations for the spin of the form >= s(s + l)h^|'0 >, 
where the right hand side is also a function of h. It turns out that when performing the limit 
h —>• 0 we get [C, D] = 0 and also >= 0, i.e., the physics of compatible observables of 
spinless systems. 

With this analysis we see that the limit h —0 of QM is in fact Newtonian Mechanics 
and not the whole domain of CM. Perhaps one of the reasons for the identihcation of CM 
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with Newtonian Mechanics in the mentioned references lies in the fact that we are used to 


work in QM with spinning systems since the early days of the quantum theory while the CM 
of spinning particles is still waiting for a complete development and improvement at least 
equivalent to the one we have achieved in the quantum domain. Even we can remember here 
that for many years, spin has been considered by physicists a strict quantum mechanical and 
relativistic property of the electron, as was pointed out by Levy-Leblond’s detailed account 
P] where the relevant references on this matter can be found. It is not strange that the 
recent history of physics had forgotten and considered unexistant the (s 7 ^ 0) region of CM. 
The spin is neither relativistic nor a quantum mechanical property of the electron. The only 
quantum mechanical aspect of the electron spin is that it is quantized. 

In QM spin is usually defined as the translationally invariant part of the total angular 
momentum operator, or in an equivalent way, as the total angular momentum measured in 
the center of mass frame of the particle, when this frame is available. In the description 
of photons we have no center of mass frame but the hrst dehnition still works. But this 
dehnition is related to the representation of the generators of the rotation group on the 
Hilbert space of states. Even more, the QM spin dynamics is not postulated, it is contained 
in its commutation relations with the total Hamiltonian. 

However, in CM spin is usually dehned by assuming new degrees of freedom, the spin 
degrees of freedom, and even a plausible dynamical equation for the spin evolution. Let us 
mention the case of the Bargmann, Michel and Telegdi equation |^. 

Following with the above scheme of Fig. we see that physics first entered into the 
QM (s 7 ^ 0) region by hand. Remember the Pauli method of introducing cr-spin matrices 
into Schroedinger’s equation, or for example when Dirac performs the linearization of Klein- 
Gordon operator to obtain the quantum mechanical relativistic description of the electron. 
But it is after Wigner’s work P on the inhomogeneous Lorentz group representations, that 
by Ending the difierent representations of the Poincare group, we can obtain the description 
of particles of any spin. But Wigner’s work not only develops the QM of spinning systems 
but it also provides a very precise mathematical definition of the concept of elementary 
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particle at the quantum level. An elementary particle is defined as a quantum mechanical 
system whose Hilbert space of pure states is the representation space of a projective unitary 
irreducible representation of the Poincare group. 

The very important expression with physical consequences, of the above mathematical 
definition lies in the term irreducible. Mathematically it means that the Hilbert space is an 
invariant vector space under the group action and that it has no other invariant subspaces. 
But it also means that there are no other states for a single particle that those that can be 
obtained by just taking any arbitrary vector state, form all its posible images in the different 
inertial frames and finally produce all linear combinations of these vectors. 

We see that starting from a single state and by a simple change of inertial observer, 
we obtain the state of the particle described in this new frame. Then let us take the 
orthogonal part of this vector to the previous one and normalize it. Repeat this operation 
with another kinematical transformation acting on the same first state followed by the 
corresponding orthonormalization procedure, as many times as neccesary to finally obtain a 
complete orthonormal basis of the whole Hilbert space of states. All states in this basis are 
characterized by the physical parameters that define the first state and a countable collection 
of group transformations of the kinematical group G. And this can be done starting from 
any arbitrary state. 

This idea allows us to define a concept of physical equivalence among states of any 
arbitrary quantum mechanical system in the following way: Two states are said to be 
physically equivalent if they can produce by the above method an orthonormal basis of 
the same Hilbert subspace, or in an equivalent way, if they belong to the same invariant 
subspace under the group action. It is easy to see that this is an equivalence relation. 
But if the representation is irreducible, all states are equivalent as basic pieces of physical 
information for describing the elementary system. There is one and only one single piece 
of basic physical information to describe an elementary object. That is what the term 
elementary means. 

But this definition of elementary particle is a pure group theoretical one. The only 


5 



quantum mechanical ingredient of it is that the group operates on a Hilbert space. Then 
one question arises. Can we translate this QM dehnition into the classical domain and then 
obtain an equivalent group theoretical dehnition for a classical elementary particle? If this 
is the case we will be able to develop the (s 7 ^ 0) region of CM in a more equivalent way. 
The answer is affirmative. 

In CM we have no vector space structure to describe the states of a system. What we 
have are manifolds of points where each point represents either the conhguration state, the 
kinematical state or the phase state of the system depending on which manifold we work. 
But the idea that any point that represents the state of an elementary particle is physically 
equivalent to any other, is in fact the very mathematical concept of homogeneity of the 
manifold under the corresponding group action. In this way, the irreducibility assumption 
of the QM dehnition is translated into the realm of CM in the concept of homogeneity of 
the corresponding manifold under the Poincare group or any other kinematical group we 
consider as the symmetry group of the theory. But, what manifold? Conhguration space? 
Phase space? The answer as has been shown in previous papers i, 0, §, is that the 
appropriate manifold is the kinematical space. 

In previous attempts, Bacry |^, considered the phase space as the homogeneous manifold 
to describe elementary objects arriving at the conclusion that the most general elementary 
particle is a system of four degrees of freedom. Three represent the position of the particle, 
being the linear momentum their conjugate variables, and the fourth is an angle a whose 
conjugate momentum is a spin component Sa such that the spin components are expressed 
in terms of these two variables and an invariant value S, the absolute value of the spin, in 
the form: 


Sx = {S^ - cos a, Sy = {S^ - sin a, 


This result was generalized independently by Kirillov, Kostant and Souriau and is known as 
the KKS theorem, showing that the coadjoint action of any Lie group dehnes on its orbits 
a symplectic structure 0 ]. 
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But the phase space, although interpreted as the state space of CM does not play the 
same role as the Hilbert space in QM at least as the dynamics is concerned. In QM the 
dynamics is stated in terms of initial \'ipi > and hnal |'0/ > states, such that the probability 
amplitude for the dynamical process \'ipi >—*• \'ipf > is given by the corresponding matrix 
element of the scattering operator < >• But both |'^j > and |V’/ > are elements of 

the same Hilbert space that at the same time it plays the role of the space that describes all 
the particle states, it also represents the kinematical space where the dynamics is running. 

In the Lagrangian approach of CM the kinematical space is the manifold spanned by 
the initial (or hnal) boundary variables that are held hxed in the corresponding variational 
formulation. It is this manifold X where the dynamics is developed, where when quantizing 
the system we obtain the natural link between the classical and quantum formalisms through 
the Feynman’s path integral approach, as shown in ||^. It is this manifold the natural base 
space to dehne on it the Hilbert space structure of the quantized system. In a formal way 
we can say that each point x ^ X, of the kinematical space of the Lagrangian formalism 
that represents the inital or final kinematical state, smears out and is transformed through 
the Feynman’s quantization into the squared integrable function with support around 
that point x and representing the particle wave function. 

Usually the Lagrangian of any classical system is restricted to depend only on the hrst 
order derivative of each of the variables g* that represent the different independent degrees 
of freedom, or equivalently, that the g* satisfy second order differential equations. But at 
this stage, if we do not know what are the basic variables we need to describe our elementary 
system how can we state that they necessarily satisfy second order differential equations? 
If some of the degrees of freedom, say gi, g2 and gs, represent the center of mass position of 
the system, Newtonian mechanics implies that in this particular case L will depend on the 
hrst order derivatives of these three variables. But what about other degrees of freedom? 
It is the homogenenity condition on the kinematical space, as the mathematical statement 
of elementarity, that will restrict the dependence of the Lagrangian on these higher order 
derivatives. It is the dehnition of elementary particle that will suply the structure of the 
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Lagrangian. 


A classical elementary particle is defined as a Lagrangian classical system whose kine- 
matical space X, is a homogeneous space of the corresponding kinematical group G. |^-§[ 

Since the Galilei and Poincare groups are 10-parameter Lie groups, the largest homoge¬ 
neous space we can hnd is a 10-dimensional manifold whose variables share the same domain 
and dimensions like the variables we use to parameterize the group. Both groups are pa¬ 
rameterized in terms of the following variables (6, a, v, a) with domains and dimensions 
respectively like h & R that represents the time parameter of the time translation, a G 
the three spacial coordinates for the space translation, v G the three components of the 
relative velocity of the inertial observers, that are restricted to u < c in the Poincare case, 
and hnally a G S' 0 ( 3 ) are three angular variables that characterize the relative orientation 
of the corresponding Cartessian frames and whose compact domain is expressed in terms of 
a suitable parametrization of the rotation group. 

In this way the maximum number of kinematical variables is also 10 , and we represent 
them hy X = {t, r, u, a) with the same domains and dimensions as above and are interpreted 
respectively as the time, position, velocity and orientation of the particle. This is the same 
description of the initial and hnal states as the relativistic spherical top in the Hanson 
and Regge canonical approach |^, but in that work different constraints are used to avoid 
dependence of the Lagrangian on second order derivatives. 

In our approach, since the Lagrangian must also depend on the next order derivatives 
of the kinematical variables we arrive at the conclusion that L must also depend on the 
acceleration and angular velocity of the particle. The particle is a system of six degrees 
of freedom, three, r, represent the position of a point and other three a, its orientation in 
space, which we can visualize by assuming a system of three orthogonal unit vectors linked 
to point r as a body frame, but the Lagrangian may depend up to the second time derivative 
of r, or acceleration of that point, and on the hrst derivative of a, i.e., the angular veloc¬ 
ity. By this dehnition it is the kinematical group G that represents the special Relativity 
Principle that completely determines the structure of the Lagrangian that represents a clas- 




sical elementary particle Point particles are particular cases of the above definition 

and their kinematical spaces are just the quotient structures between the group G and the 
homogeneous transformations subgroup of rotations and boosts, and thus their kinematical 
variables reduce only to time and position (t, r). 

We began this section by considering the classical limit of QM and we have remarked the 
different treatment the recent history of physics has done to the (s 7^ 0 ) region in both QM 
and CM formalisms. But this difference between the two approaches might have physical 
consequences. For instance, let us consider the well known QM diffraction experiment of 
a beam of particles by a certain slit. If we send a beam of point particles we expect in 
the screen a replica of the slit shape and size but if we obtain a diffraction pattern we 
interpret this outcome as a pure quantum mechanical effect. This interpretation is correct 
if for instance we perform the experiment with pions, i.e., spinless particles. But if the 
experiment is done with photons, electrons, neutrons, etc. as is usually the case, all of 
them are spinning particles. Since the spin of these particles is of order of h not the whole 
outcome of the experiment of order h is a quantum mechanical effect, because part of it 
must be related to the spin structure of the particles and must be analyzed theoretically. 
What is a quantum mechanical effect is the difference between the real outcome and that 
part that can be interpreted in classical terms but taking into account the corresponding 
classical spin contribution. A deeper understanding of the classical spin description will help 
to better understand the quantum mechanical effects. 

In this work, and following with the above argument we shall consider the classical 
contribution to the tunnel effect. Classically, a point particle conhned into a potential well 
can never escape from it if its energy is less than the top of the potential barrier. In QM there 
is a nonvanishing probability of crossing. However, if the particle has spin and this spin is of 
orbital nature as will be described in the next section, we shall see in Section V that because 
of the spin, the classical spinning particle can cross the barrier even with kinetic energies 
below the potential energy barrier, and that this contribution is not in general negligible 
when compared with the quantum tunnel effect. 
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III. CLASSICAL PARTICLE WITH (ANTI)-ORBITAL SPIN 


We are going to analyze in detail a non-relativistic Lagrangian that although it does not 
give rise to a spin 1/2 particle when quantized (see 0), however it behaves in most of the 
features simmilarly as the relativistic electron described in the mentioned reference, but the 
mathematics involved are simpler than in the general case, and illustrates the spin structure 
and dynamics for the example we want to work out in section V. 

Let Q be the Galilei group. Let us consider a Galilei particle whose kinematical space 
is A = Q/SO{?>), so that any point x E X can be characterized by the seven variables 
X = (f,r,u), u = dr/dt, which are interpreted as the time, position and velocity of the 
particle respectively. In this example we have no orientation variables and because of this 
we cannot have, when quantizing the system, half integer spin values. The Lagrangian 
will also depend on the next order derivatives, i.e., on the acceleration of the particle. 
Rotation and translation invariance implies that L will be a function of only u^, (du/dt)"^ 
and u-du/ dt = d{v? / 2 ) /df, but being this last term a total derivative it will not be considered 
here. 

Let us assume that our elementary system is represented by the following Lagrangian 

m (m /d^r\^ 


where m is the mass and the parameter u of dimensions of s“^ represents an internal fre¬ 
quency. In terms of the kinematical variables and their derivatives, and in terms of some 
group invariant evolution parameter r, the Lagrangian can be written as 

^ m f ^ m 

~ TT ~ ^T’ 


where the dot means r—derivative. If we consider that the evolution parameter is dimen¬ 
sionless, all terms in the Lagrangian have dimensions of action. Because the Lagrangian is a 
homogeneous function of hrst degree in terms of the derivatives of the kinematical variables, 
L can also be written as 
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L = —Tt + R ■ r + U ■ ii, 


( 3 ) 


where the functions accompanying the derivatives of the kinematical variables are defined 


and explicitely given by 


T 

R 

U 


dL m / dr \ ^ 

di 2 \dt) 

dL dr 

dL m d^Y 

d\i dt"^ 


m / d?r\^ 
^ \^) ’ 


(4) 


Dynamical equations are: 


1 d'^r d^r 
dt^ dt^ ’ 


(5) 


whose general solution is: 


r(t) = A + Bt + C cos cut + D sin cut, 


( 6 ) 


in terms of the 12 integration constants A, B, C and D. 

When applying Noether’s theorem to the invariance of dynamical equations under the 
Galilei group, the corresponding constants of the motion can be written in terms of the 
above functions in the form: 


Energy 
Linear Momentum 
Galilei or Static Momentum 
Angular Momentum 


H = T - 

P = R 


dU m 9 cu^ ^^9 

u ■ —— = —-U 

dt 2 2m 

dU dU 

■ —— = mu-—, 

dt dt 


u 


dt ’ 


K = mr — Pt — U, 

J = rxP + uxU = L + S. 


(7) 

( 8 ) 
(9) 

( 10 ) 


It is the presence of the U function that distinguishes the features of this system with respect 
to the point particle case. We see that the total linear momentum is not lying along the 
direction of the velocity u. The translationally invariant part of the total angular momentum 
J is the spin of the system S. It is the dependence on the acceleration that in this system 
of three degrees of freedom leads to the spin structure, such that the spin takes the form: 
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( 11 ) 


S = u X U 


m d?v dr 
a;^ dt'^ dt' 


is always orthogonal to the velocity and acceleration of point r, and its dynamics is given by 


taking the time derivative of J in (p!0D and taking into account that P is constant it leads to 


dS 

dt 


P X u. 


( 12 ) 


The spin is not a constant of the motion except in the center of mass frame, it precesses 
along the linear momentum P such that its projection S-P, called helicity, remains constant. 

If we substitute the general solution (|^) in (|7|-|l0D we see in fact that the integration 
constants are related to the above conserved quantities 


2 ^^2 , J32 


H = —B" - 
2 2 


-(C^ 


P = mB, 

K = mA, 

J = A X mB — muC x D. 


But at the same time that the particle acquires its spin we see that the 
K differs from the point particle case in the term —U, such that if we 
X = U/m, then K = 0 leads from to the equation: 

d(r — x) 


P = m- 


dt 


static momentum 
dehne the vector 


and q = r — X, dehnes the position of the center of mass of the particle that is a different 
point than r and is given by 


q = r 


u 


1 d^r 


(13) 


such that in terms of it, dynamical equations can be separated into the form: 

d^q 


df2 


= 0 , 


d^r 2/ 

^ + u.(r-q) = 0. 


(14) 

(15) 
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In terms of these variables the spin takes the form 


S = —(r — q) X mu, (16) 

and it has the formal expresion of the opposite of an orbital angular momentum with respect 
to point q of a point mass m at point r moving at the speed u. We say that this spin is of 
(anti)-orbital nature. 

Point q moves in a straight trajectory at constant velocity while the motion of point r is 
an isotropic harmonic oscillator of angular frequency uj around point q. But if q represents 
the center of mass position then, what position does point r represent? Point r represents 
the charge of the particle position. This can be seen by considering some interaction with 
an external held. The homogenenity condition of the Lagrangian in terms of the derivatives 
of the kinematical variables leads us to consider an interaction term of the form 

L/= —e0(f, r)i + eA(f, r) ■ r, (17) 

linear in the derivatives and where the external potentials are functions of t and r. We can 
also consider more general interaction terms depending on u and u, but this will not be an 
interaction with an external electromagnetic held. Then the dynamical equations are 

1 (Pr e 

where the electric E and magnetic held B are expressed in terms of the potentials in the usual 
form, and because the interaction term does not ahect the dependence of the Lagrangian on 
u, the spin and the center of mass dehnitions ( pdj) and ([T3|) respectively, remain the same 
as in the free particle case. Then dynamical equations can again be separated in the form 

^ — (E(f,r) + ux B(f,r)), (19) 

dr m 

^ + o;2(r - q) = 0, (20) 

where the center of mass q satishes Newton’s equations under the action of the total external 
Lorentz force, while point r still satishes the isotropic harmonic motion of angular frequency 
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uj around point q. This internal motion remains unaffected by the external interaction. But 
the external force and the fields are defined at point r and not at point q. Point r clearly 
represents the point charge position. In fact, this minimal coupling we have considered is the 
coupling of the electromagnetic potentials with the particle current, that in the relativistic 
case can be written as A^j^, but the current is associated to the motion of point r. This 
charge has an oscillatory motion of very high frequency u that in the case of the relativistic 
electron (see 0,||), is a; = 2mc^/h ~ 10^°s“^. The average position of the charge is the 
center of mass, but it is its internal rotational motion, usually known as the zitterbewegung, 
that gives rise to the spin structure and to the magnetic properties of the particle. 

When analyzed in the center of mass observer (see Fig. the system reduces to a point 
charge whose motion is in general an ellipse but if we choose C = D, and C ■ D = 0, it 
reduces to a circle of radius a = C = D, orthogonal to the spin that is constant in this 
frame. Then if the particle has charge e, it has a constant magnetic moment 

1 . e e 

u = -rxi = -rxu =-b, 

^ 2 2 2m 

opposite to the spin direction if e > 0, and a nonvanishing oscilating electric dipole p = er, 
orthogonal to both /i and S in this frame, such that its average value vanishes for times 
larger than the natural period of this internal motion. 

In this case the above constants of the motion are no longer conserved quantities. For 
instance, the linear momentum satisfies 


dt 


where F is the external Lorentz force. If we take the time derivative of the total angular 
momentum it leads to 


dJ 


dP dr 
dt 



dS 


( 21 ) 



The spin dynamics is given by 



dt o;^ dt^ dt 


( 22 ) 
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and remains the same expression as in the free case ([T^) , bnt now P is not a constant 


fnnction. If we snbstitnte this in the previons expression (pT]), it tnrns ont that the total 
angnlar momentnm satishes 


dJ 

dt 


r X F, 


(23) 


where on the right hand side we have the total external torqne prodnced by the helds. 
Althongh the system has magnetic moment associated to the charge motion, becanse of 
the minimal conpling interaction, the external electromagnetic held prodnces forces, and 
the corresponding torqnes with respect to the origin, bnt no other torqnes, so that a spin 
evolntion eqnation of the form 


dt 


/X X B, 


as in Ref. [Q, does not arise in this model in which the electromagnetic strnctnre of the 
particle is jnst that of a point charge withont any intrinsic magnetic moment, as it happens 
to be the case of the relativistic electron, where the charge motion in the center of mass 
frame is a circle of radins a = S/mc = h/2mc at the speed of light c. Only for anomalons 
magnetic moments we should have to consider the addition of external torques of the above 
form, where /x will be the anomalous magnetic moment. But in that case we need to have 
in the interaction Lagrangian terms of the form M^y/2, where M^i, is the dipole tensor 
that describes the intrinsic dipole structure of the particle. 

If we consider that the mechanical energy of the particle is still given by expression (^, 

then because of the interaction it changes in the form 

dH ^ ^ 

—— = u ■ F = eu • E. 

dt 


The change of the mechanical energy of the particle is the work of the external force along 
the charge trajectory. If the electric held is conservative then the change of this energy 
between two arbitrary points of the trajectory is 


H(t2) - H(ti) = e(j){ti, ri) - 60(^2, r2) 
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= e (0(ti, qi) - 0(^2, q2)) + e (0(ti, ri) - 0(ti, qi)) + e (0(^2, q2) - 0(^2, r2)) 
~ e0(ti,qi) - e0(t2,q2), 


if the difference of potential between the center of mass and center of charge position at any 
instant can be considered negligible if the external held is of smooth variation. 

Once we introduce the classical path into the free Lagrangian and integrate from the 
initial to the hnal point, we obtain the action function for the free system along this path 
that can be written in terms of the boundary kinematical variables in the form 

1 


A(fi,ri,Ui;f2,r2,U2) = 


m 


00 


-(u2 - Ui)^ sina;(f2 - 


+ a;(r2 - ri)^ sina;(f2 - ^i) 

- (h - h) (u2 - ui)^ cosoo{t2 - h) 

+ 2(^2 - h) (u2 ■ ui) (1 - cosa;(f2 - ^i)) 

- 2(r2 - ri) ■ (ui + U 2 ) (1 - cosuj(t2 - h)) 


where A = 2(1 — cosoo{t2 — ti)) — oo{t2 — fi) sina;(f2 — ^i)- This is the action function that 
must be considered when analyzing Feynman’s kernel in the path integral approach of this 
system. 


IV. CANONICAL FORMALISM 


Although the Lagrangian depends on second order derivatives we can develop the corre¬ 
sponding canonical formalism. Starting from the Lagrangian 


L = 



m 






where now the dot means time derivative, we have in this case six generalized coordinates 
that are not the independent degrees of freedom but they are the three degrees of freedom 
qi = r and their hrst derivatives q2 = dr/dt, such that the conjugate momenta are dehned 


by 0 
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_dL d (dL\ _ d\] 

dr dt\drj dt ' dr 

The phase space is a 12 -dimensional manifold and the Hamiltonian is in fact the total energy 
written in terms of the canonical variables 

2 

7 TL UJ 

^ = Pi ■ qi + P2 ■ q2 - ^ = Pi ■ q2 - — q2 - ^ P2- 


Hamilton-Jacobi equations are 

dH 

qi = ^ = q2, Pi 

( 9 pi 


dH 

dqi 


0 , 


dH dH 

q2 = =-P2, P2 = = -Pi + mq2- 

ap2 m dq2 

The ten Noether constants of the motion become in this formalism the generating func¬ 
tions of the corresponding canonical transformations of time and space translations, pure 
Galilei transformations and rotations and are explicitely given by 

TT m 2 2 

p = Pi, 

K = mqi - Pit - P 2 , 

J = qi X Pi -h q 2 X p 2 , 


and since the Poisson bracket of two constants of the motion is again a constant of the 
motion we obtain that the above constants of the motion satisfy the commutation relations 

[t/j, the] (^iklJli Pfc] ^iklPli P-k\ ^iklP-h \Jii 0, 

[Pi, Pfc] = 0 , [Pi, H] = 0 , [Ki, Kk] = 0 , [Ki, H] = Pi, [Ki, Pj] = m 5 ij, 

where [.,.] is the corresponding Poisson bracket and because [P'i,Pj] 7^ 0 they are not the 
commutations relations of the Galilei group but rather those of the extended Galilei group 
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r 3 | . Although Ki satisfies [Ki, H] = Pi, it is a constant of the motion because as we see in 
(D that K is time dependent its total time derivative is 


dK dK 

^=[K.ffl + ^ = Pi-pi = 0. 

The spin observable S = uxU = q2Xp2 satisfies the commutation relations 

Sk] = e^uSi, [Ji, = eikiSi, Pu] = 0 , 

\Si, Kj^ (^ijkUk ^ijkP2ki [*5*1 > Tf] ~(p ^ 

showing respectively that it is an angular momentum, that transforms like a vector under 
rotations, it is invariant under space translations but not under pure Galilei transformations 
and is not a constant of the motion, but satisfies the dynamical equations 0- We can 
check that the two constants of the motion 


p2 / 1 \^ 

E = H -, and = J-KxP , 

2 m \ m J 


( 24 ) 


commute with the above ten generators and they are Galilei invariant properties of the 
particle, that together with the mass m completely characterize the structure of this particle. 
They are the three functionally independent Gasimir invariants of the extended Galilei group. 
In fact E is the Galilei internal energy of the particle and Z'^ is the squared of the total 


angular momentum or spin for the center of mass observer |13 . 


V. CLASSICAL TUNNEL EFFECT 

Let us consider a spinning particle as described above under the influence of a poten¬ 
tial barrier. Sharp walls correspond classically to infinite forces so that we shall consider 
potentials that give rise to finite forces like those of the shape depicted in Fig. where Vq 
represents the top of the potential. Then the external force E{x), is constant and directed 
leftwards in the region x G (—a, 0) and directed rightwards for x G (0, b), vanishing outside 
these regions. 
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Let us assume for simplicity that the spin is pointing up or down in the 2 : direction such 
that the point charge motion takes place in the XOY plane. Let Qy and = 0, be the 
coordinates of the center of mass and x, y and z = 0, the position of the charge. 

Dynamical equations are 


d'^Qx 

dt‘^ 



d^qy 

df^ 


0 , 


(25) 


d^x 

dt"^ 


+ u‘^{x 


Qx) = 0 , 


dh 

dt"^ 


+ - in) = 0 . 


(26) 


where 


F{x) = 


—eVo/a, for x G (—a, 0), 
eVo/b, for a; e (0,6), 

0, otherwise. 


To make the corresponding numerical analysis we shall dehne different dimensionless 
variables. Let R be the average separation between the center of charge and center of mass. 
Then we define the new dimensionless position variables: 


qx = qx/R, qy = qy/R, X = x/R, y = y/R, d = a/R, b = b/R. 


The new dimensionless time variable a = out is just the phase of the internal motion, 
such that dynamical equations become 


d^qx 

da^ 


A{x), 


d^qy 

da"^ 


0 , 


d^x ^ ^ d'^y ^ ^ 


where A{x) is given by 


A{x) 


—eVo/amoj‘^R^, for x G (—a, 0), 
' eVo/bmuj‘^R^, for x G (0,6), 

0, otherwise. 
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In the case of the electron, the internal velocity of the charge is uiR = c, (see 0 and 0), 
so that the parameter ejmc^ = 1.9569 x such that for potentials of order of 100 

volts we can take the dimensionless parameter eVo/muj'^E? ~ 2 x 10“^. 

If we choose as initial conditions for the center of mass motion 

gy(0) = 0, dqy{0)/da = 0, 

then the center of mass is moving along OX axis. Then the above system reduces to the 
analysis of the one-dimensional motion where the only variables are Qx and x. Let us call 
from now on these variables q and x respectively and remove all hats from the dimensionless 
variables. Then the dynamical equations to be solved numerically are just 


d?q , , , d'^x 

= A{x), -j-^+x-q = 0, 


(27) 


da'^ da^ 

where A{x) is given by 

/ 

—2a“^10“^, for x G (—a, 0), 

26“^10“'^, for X G (0,6), 

0, otherwise. 

Numerical integration has been performed by means of the computer package ODE Work¬ 


A{x) = < 


(28) 


bench |IJ]. The quality of the numerical results is tested by using the different integration 
schemes this program allows, ranging from the very stable embedded Runge-Kutta code of 
eight order due to Dormand and Prince to very fast extrapolation routines. All codes have 
adaptive step size control and we check that smaller tolerances do not change the results. 

With a = 6 = 1, and in energy units such that the top of the barrier is 1, if we take an 
initial kinetic energy K below this threshold, K = mg(0)^/2eVo = 0.41 we obtain for the 
center of mass motion the graphic depicted in Fig. ^ where it is shown the variation of the 
kinetic energy of the particle K{q), with the center of mass position during the crossing of 
the barrier. There is always crossing with a kinetic energy above this value. In Fig. ^ the 
same graphical evolution with a = 1 and 6 = 10 and K = 0.8055 for a potential of 10^ Volts 
in which the different stages in the evolution are more evident. Below the initial values 
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for the kinetic energy of 0.4 and 0.80 respectively, the particle does not cross the potential 
barrier and it is rejected backwards. 

If in both examples the parameter a is ranged from 1 to 0.05, thus making the left slope 
sharper, there is no appreciable change in the crossing energy, so that with a = 1 held hxed 
we can compute the minimum crossing energies for different b values, Kc{b). Initial center of 
mass position has been ranged from g(0) = —3.5 to —3.0, while the center of charge initial 
position and velocity from a:(0) = —2.5 to —2.0, and x(0) = g(0), respectively. 

To compare this model with the quantum tunnel effect, let us consider a point particle 
of mass m and charge e, under the same one-dimensional potential depicted in Fig. If the 
particle of total energy E, is initially on the left hand side of the barrier, the wave function 
of this system is: 


Akx 


+ Re 


—ikx 


for X < —a. 


'0(x) = < 


C'iAi(Zi)(l — G + x/a)) + C' 2 Bi(Zi)(l — G -f- x/a)), for —a < x < 0, 
G 3 Ai(L(l — G — x/b)) -|- G 4 Bi(L(l — G — x/b)), for 0 < x < 6, 

for X > 6, 

where x is the same dimensionless position variable as before, and the constants 


(29) 


k = 


V 2mc^ ) 


y 2mc^ J 


y 2mc^ J 


eld)’ 


(30) 


and Ai(x) and Bi(x) are the Airy functions of x. The six integration constants R, T, and 
Gi, i = 1,2, 3,4, can be obtained by assuming continuity of the functions and their hrst order 
derivatives at the separation points of the different regions. The coefficient represents 
the probability of the particle to be reflected by the potential and |Tp its probability of 
crossing. 

Computed the T variable for a = 1 and different values of the potential width b, and 
for energies below the top of the barrier eVo we show on Fig. the average probability for 
quantum tunneling for four different extraction potentials Vq of 10^, 10^, 10^ and 10^ Volts. 
This average probabibility has been computed by assuming that on the left of the barrier 
there is a uniform distribution of particles of energies below eVo. 
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Let us consider for the classical spinning particle the same uniform distribution of par¬ 
ticles. Then, the function P{b) = 1 — Kc{b), where Kc{b) is the minimum kinetic energy 
for crossing computed before, represents the ratio of the electrons that with kinetic energy 
below the top of the potential, cross the barrier because of the spin contribution. This func¬ 
tion is also depicted in Fig. |^, and clearly has a decreasing behaviour with the distance b. 
The different barriers we have considered so far, can be interpreted as the one-dimensional 
extraction potentials of a probe tip placed at a variable distance b and at a constant dif¬ 
ference of potential Vq with respect to the sample. We see that for the different potentials 
shown in that hgure this probability of crossing is greater in the quantum case than in the 
classical spin contribution, but they are geting closer for higher potentials. 

The main feature of this model that contributes to tunneling is that while the center of 
mass is still on the left hand part of the barrier and since the cyclic internal motion of the 
charge is unchanged by the interaction this implies that once the charge penetrates into the 
left slope of the potential during a fraction of its cycle the kinetic energy is decreased, while 
during the remaining time of its cycle it is unaffected by external forces so that the center 
of mass motion is a kind of short decceleration periods under a constant force with short 
periods of constant velocity in between as is clearly visible in the saw-teeth shape of the 
particle kinetic energy of Fig. It turns out that the average deccelerating force acting on 
the center of mass of the particle is smaller than in the case of a pure point particle. When 
the particle has completely penetrated into the potential and not all its kinetic energy has 
been exhausted, then if the penetration is sufficient such that the charge position can reach 
during part of its cycle the right slope of the potential, it turns out that during that time the 
center of mass is under a force directed to the right and thus the kinetic energy is increased. 
It is clear that at least one of the regions (—a, 0) or (0, b) must be of the order of the radius 
of the internal motion to have an appreciable crossing. In the case of the electron model 
(see [0 and [Q) this radius is h/2mc, i.e., half the Compton wave-length of the electron, so 
that with a or 6 ot this order of magnitude as has been considered in the above examples, 
we can have a nonvanishing crossing. 
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We see that the separation between the center of mass and center of charge that gives 
rise to the spin structure of this particle model justihes that this system can cross a potential 
barrier even if its kinetic energy is below the top of the potential. 
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FIGURES 


FIG. 1. Classical limit of Quantum Mechanics. 

FIG. 2. Charge motion in the Center of Mass frame. 

FIG. 3. Potential barrier to be crossed by spinning particles. 

FIG. 4. Evolution of the kinetic energy during the crossing of the barrier with a = b = 1, initial 
energy K = 0.41 and extraction potential Vq = 100 Volts. 

FIG. 5. Evolution of the kinetic energy during the crossing of the barrier with a = 1, b = 10, 
initial energy K = 0.91 and extraction potential Vq = 1000 Volts. 

EIG. 6. Classical probability of tunneling P{b) and quantum tunneling for four different ex¬ 
traction potentials Vq- 
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